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The Mathisson equations under the Frenkel-Mathisson supplementary 
condition are studied in a Schwarzschild field. The choice of solutions, 
which describe the motions of the proper center of mass of a spinning test 
particle, is discussed, and the calculation procedure for highly relativistic 
motions is proposed. The very motions are important for astrophysics while 
investigating possible effects of the gravitational spin-orbit interaction on 
the particle's world line and trajectory. 

PACS numbers: 04.20.-q, 95.30.Sf 

1. Introduction 

70 years ago Myron Mathisson has presented the equations describing 
the motions of a spinning test particle in a gravitational field [1] 
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as as as 

where u x is the 4-velocity of a spinning particle, S^ u is the antisymmetric 
tensor of spin, m and D/ds are, respectively, the mass and the covariant 
derivative with respect to the proper time s; R^ pa is the Riemann curvature 
tensor of the spacetime. (Throughout this paper we use units c = G = 1. 
Greek indices run 1, 2, 3, 4 and Latin indices 1, 2, 3; the signature of the 
metric (-,-,-,+) is chosen.) Equations (P), ([2]) were supplemented by the 
condition [1] 

S^Uu = 0. (3) 
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(In special relativity condition ([3]) was introduced by J. Frenkel |2j.) It is 
known [31 Hj that in the Minkowski spacetime the Mathisson equations ([!])- 
([3]) have, in addition to usual solutions describing the straight worldlines, a 
family of solutions describing the oscillatory (helical) worldlines (as a partial 
case, this family contains the circular solutions.) 

In [3] the oscillatory solutions of equations ([I])-© were connected with 
" Zitterbewegung" . C. Moller proposed another interpretation of these so- 
lutions: he pointed out that 1) in relativity the position of the center of 
mass of a rotating body depends on the frame of reference, and 2) condition 
([3]) is common for the so-called proper and non-proper centers of mass [5j. 
The usual solutions describe the motion of the proper center of mass of a 
spinning body (particle), and the helical solutions describe the motions of 
the family of the non-proper centers of mass [5] . 

Later A. Papapetrou derived equations (H|), © by the method which 
differs from Mathisson's one [B], and, instead of ([3D, the non-covariant con- 
dition 

S l4 = (4) 

was used [7] in concrete calculations. 

To avoid the superfluous solutions of equations ([T]), ([2]), W. Tulczyjew 
and W. Dixon introduced the covariant condition [8] 

S^P V = 0, (5) 

where 

DS uft 

P" = mu v + u M ^— (6) 
as 

is the particle's 4-momentum. In contrast to relation (j3J) the Tulczyjew- 
Dixon condition ([5|) picks out the unique worldline of a spinning test particle 
in the gravitational field. That is, equations (PJ, (J2J) under ([5]) do not admit 
the oscillatory solutions. However, the question arises: is this worldline 
close, in the certain sense, to the usual (non-helical) worldline of equations 
(PQ), ((2J) under condition ([5]), for example, in a Schwarzschild field? It is 
simple to answer this question when the relation 



m\u v \ 



u, 



ds 



(7) 



takes place, because in this case condition ([5]) practically coincides with ©. 
For example, these conditions are equivalent for post-Newtonian expansions 
[9]. However, a priori another situation is possible for the highly relativistic 
spinning particle. Naturally, this situation must be investigated carefully. 

Note that very condition ([3]) was derived in some papers by different 
methods |10t \TT\ [12] , and we agree with the conclusion that this condition 
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"... arises in a natural fashion in the course of the derivation", [TT], p. 112. 
That is, condition ([3]) is necessary, though often, with high accuracy, it can 
be substituted by condition ([5|). 

We stress: the existence of the superfluous (oscillatory) solutions of 
equations ([I]), ([2]) under condition Q is not a reason to ignore this condition. 
The point of importance is that just among all solutions of equations (P)-® 
the single solution describing the motion of the particle's proper center of 
mass can be found. Obviously, it is necessary to know how this solution can 
be identified among others. 

In the focus of this paper is just the initial Mathisson equations ([!])- 
([3]) . Our purpose is to investigate non-oscillatory solutions describing highly 
relativistic motions of the spinning particle in a Schwarzschild field. 

Note that the information on all possible types of motions of the spinning 
test particles in the gravitational fields is important for astrophysics, for 
more fine investigations of the gravitational collapse and other astrophysical 
phenomena. 

This paper is organized as follows. The known integrals of the strict 
Mathisson equations (HI)-© in a Schwarzschild field, the energy and angular 
momentum, are used in section 2 for reducing the order of differentiation in 
these equations. The possible procedure of finding the values of the energy 
and momentum parameters, which correspond just to the motions of the 
proper center of mass, is discussed in section 3. This procedure is realized 
in section 4 for the motions close to highly relativistic equatorial circular 
orbits in a Schwarzschild field. In section 5 the illustrations of computer 
calculation are presented. We conclude in section 6. 

2. Mathisson equations for equatorial motions 
in a Schwarzschild field 

Let us consider equations (HI)-© for Schwarzschild 's metric in the stan- 
dard coordinates x 1 = r, x 2 = 9, x 3 = tp, x 4 = t for the equatorial 
motions of a spinning particle with spin orthogonal to the motion plane 
9 = 7r/2. Then the non-zero components of the metric tensor are: 



where M is the Schwarzschild mass. Due to the symmetry of Schwarzschild 's 
metric equations (pQ) , (|2]) have the integrals of the energy E and the angular 
momentum L which for the equatorial motions can be written as [13J 




(8) 
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L = —mU3 — <?33U^ — -933,ii • (9) 

For the equatorial motions with spin orthogonal to the motion plane 
9 = 7r/2 equations (|2|) can be solved separately from ([1]). Indeed, taking 
into account relations ([3|) and ([8]) it is not difficult to obtain from ([2]) all 
non-zero components S^ u : 

gl3 _ _Cj31 _ u 4:So !4 _ _g41 _ u sSq 

r ' r ' 

5 34 = -5 43 = -^, (10) 
r 

where So is the known constant of spin [10] 

5*0 = -S^S^. (11) 

(We stress that expressions (fit)]) satisfy all equations of set fl2J).) 

Now we shall consider equation ([1]). It is known that these equations 
under condition ([3]) contain the second order derivatives with respect to s. 
However, in our case of Schwarzschild's metric, due to the integral E and L, 
it is possible to reduce the order of differentiation by the standard procedure 
of the differential equations theory. Using (j8lH(flC)]) and the relation u^u 11 = 
1 after direct calculations (which are quite simple but rather lengthly) we 
get from equation (pQ) the two non-trivial equations for r{s) and (p(s): 

r 2 „ (. 3M\ .o rE 



+ 2r [ 1 - — ) 2 - — 
r V r J S 



1 / 3M\ L 
+ - 1 + — 

r V t J rbo 



1/2 



roi / 3M\ tp m + L(p 

Cp = + r[ 1 + — — 

r \ r J r rb^r 



r 2 + (l- — )(l +r 2 ^ 



(12) 

1/2 



(13) 

where a dot denote the usual derivatives with respect to s. (Two other 
equations of set ([I]) are satisfied identically.) So, equations (fl~2j) . (fl~3|) do not 
contain the third coordinate derivatives. However, in these equations the 
quantities E and L are present as the parameters which are not determined 
by the initial values of r, tp, r = u , (p = n 3 only. (According to ([9]), for 
the determination of E and L the second coordinate derivatives must be 
given as well.) That is, equations (|12|) . (|13p. as well as the initial Mathisson 
equations ([T]), describe the motions both of the proper center of mass and 
the non-proper centers. 
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In other words, equations (|12p . (|13|) contain the non-oscillatory and 
oscillatory solutions. The question of importance is: which values E and L 
correspond just to the proper center of mass at the arbitrary initial values 
of r, if, r, cp for a spinning test particle? It is easy to answer this question 
when the motion of such a particle is close to the geodesic motion: then 
approximately E = mu^ and L = —mu^, i.e., we write the relations for 
the geodesic motion. However, we have not any proof that the worldline of 
a spinning particle is close to the corresponding geodesic worldline for all 
physically admitted initial values of r, ip, r, (p. On the contrary, our results of 
investigations of the gravitational spin-orbit interaction in a Schwarzschild 
field [144 f!5l [T6] show that just highly relativistic motions of a spinning 
particle must be studied carefully. Therefore, in the next section we shall 
consider the procedure of choosing the values E and L for the proper center 
of mass for highly relativistic motions. 

For further calculations it is convenient to write equations (|12p . (fT3' 
terms of the non-dimensional quantities 



m 



_ s dr dip r 

T= M' y =d? Z = M d? P =M' £ ^mM' 



(14) 



(In the following we shall put e > 0, without any loss in generality.) Then 
according to equations (I12h , (I13p we have the set of the first-order differential 
equations 



dY Y 2 ( 3 
— = — + p 1-- 
dr p \ P 



2Z 2 + 



1 



fiZp 



+ 



dZ_ 



Y- \ { I --) {l + Z 2 p 2 ) 



1/2 



(15) 



YZ Z 2 + l p 2 3Z 

+ P v \ Z ^ 

P Y \ p 



where 



Y 2 + ( 1 - - ) (1 + Z 2 p 2 " 



dp 
Ih 

dcp 

ME 
"Sb"' 



Y, 



Z. 



So' 



1/2 



(16) 
(17) 
(18) 

(19) 
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For correct studying the physical conclusions of equations (I15p - (ll8p . it is 
necessary to take into account the Wald condition for a test spinning particle 
[17j . therefore we put 

£ s J? « »■ < 20 > 



Because in the following we shall compare some solutions of equations (|15D 
(|18|) with solutions of the geodesic equations in Schwarzschild's metric, it 
is useful to write here the last equations for the equatorial motions in the 
coordinates r, tp 

r J r 2 



— )--2> (21) 



2 , N 

(p = — ftp. (22) 
r 

Using notation (114j) we rewrite equations (I21j) . (j22|) as the set of the first- 
order differential equations 

— = -2—, (24) 
ar p 

plus two equations which coincide with (fT7|) , ([THD • Equations (fT5j) , (fTo^) and 
()23|) . (I24p correspondingly are essentially different. As well as equations 
(fl~2j) . (fl~3|) . the first two equations of set (fl~5j) - (fl~8|) contain the constants of 
the energy and angular momentum of a spinning particle. According to 
discussion above, different values of the parameters \i and v in (fT5|) . (fl~6j) at 
the fixed values of Y, Z, p correspond to different centers of mass, namely to 
the single proper center of mass and to the set of non-proper centers. 

3. Possible way of finding the parameters /i and v for the proper 

center of mass 

The right-hand sides of equations (fl~5|) . (fl~6j) are too complicated for 
investigations in general case. However, a priori we cannot exclude the 
possibility to show the essential difference between the motions of the proper 
and non-proper centers of mass during the short time interval after the 
beginning of the corresponding motions, namely, when the displacement of 
the values Y, Z, p, tp from their initial values Yq,Zq,pq, tpo are considered in 
the linear approximations in the quantities 

£l = — , 6 = — 7} — , 6 = • (25) 

Z p 
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(In (I25p we do not write the displacement <p — (fo because the right-hand 
sides of equations (fT5]) - (fT8j) do not depend on ip and equation (fT8j) is trivial, 
i.e., according to (|18p the value ip is determined by simple integration of Z.) 
Let us check this possibility, i.e., consider equations (fT5]) - (fT7|) in the linear 
in £ approximation. Then by direct calculations we obtain 

= (oio + a>\\v + a\ 2 p) 6 + (020 + a 2 iv + a 22 p) £2 
+ (030 + a-ii" + 032^.) £3 + aoo + a iu + a 02 p, (26) 
^P- = (ho + h\v + b 12 p- + 6 + (ho + Z>2if + &22M + &23~1 £2 



7 I " ±\J 1 "11" ' ~ _L .£ f-~ 1 "_LO I "si 1 l"^U 1 "^±" ' ** 1 ~ Z,tJ 

ar \ e / V £ . 

+ ^30 + hi" + &32/i + &33~) Cl + fo 00 + hi" + &02/U + £>03^, (27) 

-7^ = C106 + coo, (28) 
ar 

where the coefficients a, 6, c with the corresponding indexes are expressed 
through Yq,Zq,pq. We shall use the expressions of these coefficients in the 
approximation 

Z 2 pg»l, ^ 2 «^oVo- (29) 

According to notation (|14p. relations f)29[) mean that the tangential compo- 
nent of the particle's initial velocity is highly relativistic, and, in addition, 
that the tangential component is much greater than the radial component. 
First, it corresponds with our aim to investigate just highly relativistic mo- 
tions. Second, because the deviation of a spinning particle is caused by 
the gravitational spin-orbit interaction, most clearly this deviation can be 
shown in the case when the tangential velocity is dominant (in particular, 
for the circular or closer to the circular orbits [16j). 

According to the theory of differential equations, the general solution 
of linear equations (|26p - (|28p is determined by the combination of e XiT 
(i=l,2,3), where Aj are the solutions of the third-order algebraic equation 

A 3 + C 2 A 2 + CiA + C = 0. (30) 

Here the coefficients Cj (j = 0,1,2) can be expressed through a,b,c and 
depend both on po, Yq, Zq,e and on the parameters p,, v. Our task is to find 
such concrete values p, v which at the fixed po, Yq, Zq, e determine just the 
motion of the proper center of mass. We begin by analyzing the expressions 
Cj for the simple case when the worldline of a spinning particle is close to 
the corresponding geodesic worldline. 
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3.1. Expressions Cj for quasi- geodesic motions 

If the parameter e in (I20p is sufficiently small, for any fixed values 
Po,Yq,Zq the motion of the proper center of mass is close to the geodesic 
motion. Then we can write approximately E = mu^, L = —mu^ (see 
expressions ©), and according to (fl"ij) . (fT9j) we have 

1 / 2\ 1 / 2 In 
H=-\Zq\po\1 , u = -p Z . (31) 

£ \ PO/ ^ 

Let us consider expressions Cj at p, v from (|3ip . It is easy to show that all 
these expressions have a common feature: all greatest terms with the large 
value 1/e are cancelled. It means that the corresponding largest terms with 
1/e are absent in the expressions Xj as well. So, just the values p,u for 
the proper center of mass give the minimum values of \j. We stress that 
similar situation takes place for motions of the proper center of mass in the 
Minkowski spacetime: for the proper center Xj = (it corresponds to the 
straightforward motions), and for the non-proper centers Xj are proportional 
to M/Sq (the oscillatory motions). 



3.2. Expressions Cj for the highly relativistic circular motion with r = 3M 

In [16] we have considered the case of the circular motion of the proper 
center of mass in a Schwarzschild field with r = 3M. By the notation (|14p 
in this case we write 

3-3/4 

p = 3, Y = 0, Z = — . (32) 

v £ 

For values (fT9l) . from ([9]) we obtain 

^ = 3 ~l/4 £ -3/2 5 l/ = _ 3 5/4 e -3/2 > (33) 

Let us estimate the values Cj. Taking into account ([33j) it is easy to check 
that, as well as in the previous case, in all expressions Cj the largest terms 
with 1/e are cancelled. 

It is naturally to suppose that the similar feature takes place not only in 
the two partial cases above. Therefore, below we shell check do the criterion 
of excluding the largest terms with 1/e in the expressions Cj can be used 
for finding the values p and v which pick out just the motion of the proper 
center of mass. 
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4. Expressions /i. and v for the motions of the proper center of 
mass at 2 < po < 3 under condition (1291) 

As we pointed out in [16] equations CO)-© admit in a Schwarzschild 
field the solutions describing the circular equatorial orbits in the region 
2M < r < 3M. By notation (I14p for these orbits we have 



2 < p < 3, Y = 0, Z 



2 x 1/4 



P 



-1/2 



-1/2 



(34) 



The expression Z from (|34p is valid beyond the small neighborhood of the 
value p = 3. (This neighborhood was considered in |16j.) At expressions 
(I3l) . it follows from © that 



fJL ~ £ 



-1/2 



-1/2 



(35) 



It is interesting to consider the non-circular orbits which deviate from (pP 
due to lo 7^ under condition (|29|). Let us consider the expressions Cj in 
the case 

2<p <3, l«y 2 «Z 2 ^, Z = ke- l ' 2 : (36) 
where according to 



k = 1 



po 



P0 



1/4 



/50 



-1/2 



(37) 



(38) 



and fci, £)2 are some coefficients which do not depend on e. Our task is to 
find such values k\ , hi which ensure excluding the largest terms with 1/e in 
the expressions Cj. 

Using /x, v from (|37p we obtain the conditions under which the largest 
terms with 1/e in the expressions Cj are cancelled: 



2 \ l l 2 

3k 2 ( Po -3) -2 Po kk 1 -2k [ 1 £> 



Po/ 



2 \ 1 



- A; 2 (/) - 3) + p^kkx + /cfc 2 ( 1 



2 x 1/2 



k po — pokk\ — kk2 



Po 



Po 

2 \ -V2 



+ 1 



Po 



1 

Po 



1V /2 

Po) 

2_\-i/2 
Po/ 



0. 



(39) 

(40) 
(41) 
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It is easy to check that at the value k from (j37H among three linear in ki, 
hi algebraic equations ([39|) - ([^T]) there are only two independent equations. 
Then the solution of ([3"9l-(HTD is: 



1 



i 1 -- 


r 


1 - 


3 


\ Po, 






Po 


PoJ 


-1/4 


1 - 


3 




Po 



-3/2 



-3/2 



3 3 
1-- + - 

Po Po 



9 15 
1-- + — 

Po Po 



(42) 



(43) 



Relations (j4~2j) . (|43p can be used in computer integration of equations (fl"5j) 
(fl"8l) under conditions (1361). 



5. Examples of computer integration of equations (j!5|) — (j!8|) 

Graphs of p(r) and p(</?) according to equations (fl"5l) - (fl"8j) under relations 
(I42|) . (|43p are shown by the thick lines in figures 1 and 2 correspondingly. 
For comparison, the thin lines in the same figures show graphs of p(r) 
and p(ip) according to geodesic equations (f23"j) - (|2"4"j) . If relations (|4"2"j) . (|4"3"]) 
are violated, computer integration shows the oscillatory solutions. Similar 
situation takes place for orbits with po = 3. An example of the oscillatory 
solution is presented in figure 3. 



Fig. 1. Graph of p(r) according to equations p5|) - lfT8|) under relations (|42l . 
(HP at e = 10~ 6 , po = 2.5, F = 0.3 and Z determined by J37]) (thick line). 
Graph of p(r) according to geodesic equations (|2"3"|) - (f!M)) at the same p , Yq, Z 
(thin line). The horizontal line p = 2 corresponds to the horizon surface. 



Fig. 2. Graph of p(ip) in the polar coordinates according to equations (|15p - 
(fT5|) (thick line) and P3"|) ~ (f!M)) (thin line). The values po, Yq, Zq are the same as 
in Fig. 1. The circle p = 2 corresponds to the horizon surface. 



Fig. 3. An example of the oscillatory solution for p(r) according to equations 
(Unj-dlHI) at e = 10" 4 , po = 3, Y = 2.5 x 10" 3 and Z determined by (J3TJ> when 
relations ([3^)1 are slightly violated. 



6. Summary 

The computer integration of equations (fl"5"|) - (fl"8j) at relations (f4*2|) . (j33]) 
shows that these relations are suitable for choosing solutions, which describe 
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the motions of the proper center of mass. As a result, according to figures 
1 and 2, we conclude that (under those conditions indicated in the captions 
at these figures) the force of the interaction of the particle's spin with the 
gravitational field acts as the repulsive one. Due to this force the spinning 
particle falls on the horizon surface during longer time as compared to the 
corresponding particle without spin (figure 1). Moreover, according to figure 
2, considerable space separation of the corresponding spinning and non- 
spinning particles takes place within a short time, i.e., within the time of 
the particle's fall on the Schwarzschild horizon. 

There is also a possibility of generalizing the above onto case of a Kerr 
field. 

R.P. thanks Professors A. Trautman and E. Malec for useful discussions 
and hospitality in Warsaw and Cracow. 
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